Role of fluctuations in a snug-fit mechanism of KcsA channel selectivity 
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I. INTRODUCTION 

The KcsA potassium channel belongs to a class 
of K+ channels that is selective for K+ over Na+ 
at rates of K+ transport approaching the diffu- 
sion limit [Ij. This selectivity is explained ther- 
modynamically in terms of favorable partitioning 
of K+ relative to Na+ in a narrow selectivity fil- 
ter in the channel. One mechanism for selectivity 
based on the atomic structure of the KcsA chan- 
nel [2] invokes the size difference between K+ and 
Na+, and the molecular complementarity of the 
selectivity filter with the larger K+ ion An 
alternative view holds that size-based selectivity 
is precluded because atomic structural fiuctua- 
tions are greater than the size difference between 
these two ions [4\. We examine these hypotheses 
by calculating the distribution of binding ener- 
gies for Na+ and K+ in a simplified model of the 
selectivity filter of the KcsA channel. We find 
that Na+ binds strongly to the selectivity filter 
with a mean binding energy substantially lower 
than that for K+. The difference is comparable 
to the difference in hydration free energies of Na+ 
and K+ in bulk aqueous solution. Thus, the av- 
erage filter binding energies do not discriminate 
Na+ from K+ when measured from the baseline of 
the difference in bulk hydration free energies. In- 
stead, Na+/K+ discrimination can be attributed 
to scarcity of good binding configurations for Na+ 
compared to K+. That relative scarcity is quanti- 
fied as enhanced binding energy fluctuations, and 
is consistent with predicted relative constriction 
of the filter by Na+. 

A tetrameric constellation of four TTVGYG amino 
acid sequences comprises the selectivity filter shown in 
Fig. n This K+ signature sequence is strongly conserved 
across a variety of K+ channels 0. In each binding site, 
for example the site defined by (VG)4, the ion is coor- 
dinated by four carbonyl oxygens from above and four 
from below. This binding site provides a snug fit for 
K+. Thus it is argued that when K+ enters the binding 
site its dehydration is compensated by favorable interac- 
tions with the carbonyl oxygens contrast, the 
smaller Na^ ion is unable to interact optimally with all 
of the available carbonyls. This imperfect compensation 
of dehydration is suggested as the basis for the observed 
selectivity. Implicit in this hypothesis is the idea that the 



channel is stiff, and poor coordination of all the carbonyl 
oxygens leads to weaker binding of the ion to the binding 
site. We will refer to this as the snug-fit mechanism. 

A focused attempt to evaluate the snug-fit mechanism 
at a molecular level is the recent work of Noskov et al. 
iJj. They argued that size-based selectivity is precluded 
because atomic structural fluctuations are greater than 
the difference in Pauling radii for the two ions. Through 
numerical experimentation involving extensive free en- 
ergy calculations, they concluded that local interactions 
leading to structural flexibility of the binding site pro- 
vided a key to selectivity for K+. They pointed-out that 
the selectivity of the binding site is sensitive to carbonyl- 
carbonyl repulsive interactions. Artificially turning-off 
carbonyl-carbonyl electrostatic repulsion while retaining 
carbonyl-ion electrostatic attractions was found to shift 
the thermodynamic selection in favor of Na+ . 

In the analysis below we follow Noskov et al. 4] in 
studying a simplifled model of the selectivity filter. A dis- 
tinction of the present work with that of Noskov et al. is 
the use of more parsimonious statistical thermodynamic 
analyses. This sharpens the physical points that may 
be observational. On this basis we might expect experi- 
mental validation of the key result of these analyses, i.e., 
that constriction of the filter and enhanced fluctuations 
accompany each other in the binding of an under-sized 
ion such as Na^. Further experimental investigation of 
the joint association of these factors might help to clarify 
the combined roles of fluctuations and a snug flt 0, 01 in 
the mechanism of ion selectivity by the KcsA channel. 



II. THEORY 

The equilibrium selectivity of the filter can be char- 
acterized by the difference in the interaction free energy 
for transferring a Na^ ion from water into the selectivity 
filter compared to the case for a K+ ion. Thus we study 
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[/^Na+ (aq) - /^K+ (aq)] 
= A/x°'= (filter) - A^l'' 



(aq) . (1) 



Here /i^(aq) (X=K+, Na+) is the hydration free energy 
of the ion, and /i^ (filter) is the analogous quantity in the 
selectivity filter, /i^(aq), the excess chemical potential, 
is that part of the chemical potential that would vanish if 
intermolecular interactions were to be neglected. Thus, 
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fi^{aq) is understood to be referenced to the ideal gas 
result at the same density and temperature. This tech- 
nical point deserves emphasis because we could adopt a 
standard state in which Af/'^ (aq) would vanish, but that 
would not change any physical consideration. 
The potential distribution theorem 



(2) 



tells us how /i^ may be calculated. Here e is the binding 
energy of the X ion to the medium; Px (e) is the proba- 
bility density function of this interaction energy. Px (s) 
is generated with the ion and the medium fully coupled 
at temperature T = l/k^f3 where fce is Boltzmann's con- 
stant. In the present approach ion positions contributing 
to the sample are those corresponding to its natural mo- 
tion in the ion-protein system. If Px(£) is well-described 

by a gaussian of mean (e) and variance cr^ = ^(e — {e))^ 

then Px(e) cx e-(^-<'^))'/2'^' ^ and 



(3) 



The width parameter a gauges a fluctuation contribution. 

Note that in sampling the fully-coupled system the 
term /3o'^/2 raises the chemical potential above the mean 
binding energy (e). This should be contrasted with the 
case of sampling from the uncoupled subsystems. In the 
uncoupled case, contributions beyond a mean field term 
lower the free energy. The distinction reflects the fact 
that the mean binding energies are computed from differ- 
ent probability distributions. P^^ (e) is the probability 
distribution function for the binding energy when the ion 
and the medium are uncoupled, and is given by 



(4) 



Intrinsic fluctuations of binding energies are associated 
with these probability distributions, and contribute to 
the free energies in a natural way. 



{e) = j ePx {e)de 

= fi'^-kBT J Px(£)ln 



de 



(5) 



which uses Eq. 0] Here the additional fluctuation con- 
tribution, the right-most term in Eq. [SJ suggests an en- 
tropic contribution to the free energy of binding beyond 
the mean interaction energy for the joint system. Note 
specifically, however, that (e) is expected to be temper- 
ature dependent so that this additional fluctuation con- 
tribution is not an identification of the thermodynamic 
entropy contribution. 

Additional physical perspective can be obtained from 
the formal relation 



£ e 

kT In J Px{e)e^''de-kTln J P^ {e)de , (6) 



which is true independently of the binding energy cutoff 
parameter e; see for example 0, 0|. We expect Px (e) 
to be concentrated near e « (e), and thus it is natural to 
choose £ so that 



(s) = kT In J Px(e)e'3"de . 
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Then 
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The natural estimate is e 

To the extent that the gaussian estimate Eq.^jis accu- 
rate, the observed variance of binding energies observed 
for the fully coupled system teaches us about available 
states of the filter alone. Fig.|21illustrates this connection 
between the available states for the uncoupled selectivity 
filter and the observed fiuctuations for the fully coupled 
system. 

We anticipate results that follow by noting that /3cr^/2 
takes values ranging from 12 kcal/mol to 6 kcal/mol, 
roughly, for the cases of Na+ to K+. In fact, e will be 
substantially lower for Na+ than for K+. The physical 
interpretation from Eq. |H1 is that we would need sample 
sizes as big as 
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/ PPie)de 



exp[/3(M?f -(£))] 



5 X 10^ (Na+), 
2 X 10-*, (K+) 



(9) 



to find a filter configuration that provides a binding en- 
ergy e < e by sampling from uncoupled systems. If we 
found one favorable configuration in such samples, the 
probabilities would be all correct. Configurations of the 
uncoupled filter that would be favorable for Na+ are thus 
less probable than those that would be favorable for K+ . 

As a check for the gaussian approximation Eq. El we 
have also computed A/i"'' (filter) by transforming K+ to 
Na+ on the basis of a coupling-parameter integration 
through 20 intermediate states. This is, of course, a sim- 
ple algorithmic approach to evaluation of Eq. |21 Note 
that by introducing multiple gaussians |0, the single 
gaussian approximation can be refined to achieve quan- 
titative agreement with coupling-parameter integration. 
We do not pursue this point further here in favor of phys- 
ical clarity. 



III. CALCULATIONS 

The atomic coordinates for the protein structure were 
obtained from the Protein Data Bank (PDB ID: 1BL8 
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). Only residues 63 to 85 were retained in the sim- 
ulation model (Fignj. Residues 74 to 79 comprise the 
selectivity filter. Following earlier notation j^, residues 
76 (V) and 77 (G) comprise the binding site denoted as 
S2. All our calculations pertain to a single ion located in 
the S2 site. 

The N-terminal group of each of the four protein chains 
was acetylated, the C-tcrminal group was amidated, and 
all eight carboxylates were ionized. Hydrogen atom po- 
sitions were built and the structure was energy mini- 
mized keeping only the non-crystallographically deter- 
mined atomic positions free. This initial structure pro- 
vided the starting point for further simulations. 

The heavy atoms outside the selectivity filter ex- 
perience a harmonic mean field with force constant 
10 kcal/mole-A^, and similarly the heavy atoms of the 
selectivity filter felt harmonic external forces correspond- 
ing to k =10, 5, 2.5 kcal/mole-A^. We also considered a 
hybrid case in which k — 0.0 for the carbonyl oxygens of 
the S2, but k = 2.5 for all other atoms of the filter; this 
is plotted as fc = 0.0 in Fig. El 

Molecular dynamics studies were carried out with the 
NAMD program [13 using the CHARMM27 13] force- 
field. A temperature of 298 K was maintained by velocity 
scaling. The Lennard-Jones parameters for the ions were 
from hj. All the non-bonded interactions were switched 
off from 17 A to 20 A. A shorter cutoff does not affect 
the results. 

In using Eqs.ElandOl we could consider the ion-protein 
interaction in full. But the binding energy e is particu- 
larly sensitive to local, near-neighbor interactions, and 
far-field contributions are expected to be less discrimi- 
nating between Na"*" and K+ in the filter. Therefore, 
to obtain differences A/j,™ (filter) we compute only the 
local interactions between the ions and ligands explic- 
itly present. We investigated several cases for includ- 
ing ion-filter interactions of different types in evaluating 
the binding energy; our quantitative results changed by 
about 20% in the extremes, and those different approxi- 
mations did not affect our conclusions below. For clarity 
then, we present results for the case where interactions 
contributing to the binding energy are limited to the car- 
bonyl groups alone. Complementary coupling-parameter 
integrations with the same scheme for considering inter- 
actions were also performed. 

For bulk hydration studies, the SPC/E water 
model was used as much of our preliminary work had 
been done with this water model. The ion parameters 
had been developed for the TIP3P water model, but the 
hydration results with SPC/E are not significantly dif- 
ferent. The experimental partial molar volume of the ion 
|16| was used to fix the simulation volume. For aque- 
ous simulations, long range electrostatics were treated 
using Ewald summation, with non-electrostatic interac- 
tions cutoff at 8.8 A. 

The coupling-parameter transformation of K+ to Na+ 
was carried out in 20 steps, with 5 ps (10 ps for aqueous 
runs) for equilibration and 20 ps for statistical averag- 



ing. For calculations using the gaussian model Eq. |21 a 
simulation length of 0.5 ns was used and the data stored 
every 20 steps (25 steps for aqueous runs) for analysis. 

IV. RESULTS AND DISCUSSION 

A. Bulk hydration 

The bulk hydration free energies set the baseline for 
filter selectivity as is evident in Eq. ^ The coupling- 
parameter integration yields A/i™ (aq) = —20.7, in 
good agreement with The gaussian model gives 

—23.0 kcal/mole. This result is reasonable for a sin- 
gle gaussian description of aqueous hydration and 
is consistent with the coupling-parameter integration re- 
sult. 

From Eq. ^ the relative selectivity is Afi°^ — 
A^™ (filter)-|-20.7 kcal/mole. It is estimated that A/z™ is 
about 6 kcal/mole which is about three times smaller 
than the difference in bulk hydration free energies. Thus 
the hydration thermodynamic properties in bulk aque- 
ous solution alone play a dominating role in the filter 
selectivity. 

B. Filter results 

Fig.Olshows the ion-carbonyl-oxygen pair distribution 
function. Notice that carbonyl oxygen atoms approach 
the Na+ more closely, consistent with the smaller size of 
Na"*" and the more favorable coulombic interactions ob- 
tained then. The carbonyl oxygens are also more delocal- 
ized in the case of Na"*" than that of K+. These observa- 
tions imply that the selectivity filter has conformational 
fiexibility to readily accommodate Na+, as opposed to 
the static structural picture discussed in standard texts, 
for example Fig 13.25 in ^ and Fig 11.24 in Ts^, in 
which the ion moves about within a fixed filter struc- 
ture. Instead these results suggest that the Na+ would 
achieve favorable binding energies by shifting the ensem- 
ble of configurations for the filter away from that for the 
case of K+. 

This relative constriction can be directly observed as 
Fig. ^ shows. This narrowing is expected to be sensitive 
to the repulsive interactions between carbonyl-oxygen 
atom pairs which have been implicated in selectivity jj]. 
But selective manipulation of carbonyl 00 repulsions is 
artificial and, thus, experimental investigation of the ne- 
cessity of this point docs not seem likely. 

Fig. |S1 shows the distribution of ion-carbonyl binding 
energies for the k = case. It is clear that the mean in- 
teraction energy of Na"*" with the filter is much lower than 
for K+ . In addition, the distribution of this interaction 
energy for Na"*" is broader than that for K"*", consistent 
with relative widths of the distributions of Fig. O 

Fig. El collects the results for the various cases consid- 
ered. These results are not importantly sensitive to fc. 
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Consistent with experiments, A/i™ > 0: the filter selects 
K"*" over Na"*". The magnitude of Ajj,'^^ is also roughly 
consistent with the estimate of about 5 — 6 kcal/mole 

Nevertheless, Na+ achieves the energetically more- 
favorable binding to the filter. A(e), the difference 
in mean binding energy of Na+ to the filter relative 
to that for K+, is between between —22 kcal/mole 
and —25 kcal/mole, consistent with the energy differ- 
ence noted in Q. These values are about the same 
as the difference in aqueous hydration free energy of 
—23 kcal/mole. Thus, filter discrimination against Na+ 
on the basis of binding energies alone is implausible. 

However, the difference in fluctuation contributions, 
A(/3cr^/2), must also be considered, and the magnitude 
of this difference is comparable to the net selectivity, as 
Fig. shows. A(/3cr^/2) is roughly 5 — 6 kcal/mole, in- 
dicating that filter configurations conducive to binding 
Na+ are comparatively rare, consistent with Eq. |^ It is 
this difference in fiuctuation contributions for Na+ rela- 
tive to K+ that shifts the balance in favor of K+. 



V. CONCLUSIONS 

Na"*" binds strongly to the selectivity filter with a mean 
binding energy substantially lower than for K+. The 
difference is comparable to the difference between the 
hydration free energies for these two ions in bulk aque- 
ous solution. Since the ion-sorting ability of the KcsA 
K+ filter must be considered from the baseline of this 
substantial difference in bulk hydration free energies, we 
conclude that the average filter binding energies alone do 
not provide significant discrimination of Na+ from K+. 

Strong binding of the smaller Na+ also constricts the 
selectivity filter. From the point of view of the filter, dis- 
crimination against Na+ results from the fact that the ob- 
served Na~''-narrowed conformations are rare occurrences 
in the ensemble of conformations favorable to K+. This 
effect is described in the thermodynamics as an observed 
fiuctuation contribution that is destabilizing for Na^ rel- 
ative to K+. The key result here is the association in the 



case of Na+ of strong binding, constriction of the filter, 
and enhanced fluctuations. 

With respect to the snug-fit view of channel selectivity, 
it is clear that a favorable binding energy is important. 
In the case of the smaller ion (Na+), this is achieved at 
the expense of stronger energetic and positional (Fig. ^ 
fluctuations. For the ion considered to have the better 
geometrical flt (K+), this is achieved at a lower cost in 
energetic and positional fluctuations. In this sense, the 
size difference between Na^ and K"*", and the molecular 
complementarity of the selectivity filter with K+ do play 
a role in channel selectivity. From the fluctuations point 
of the view, the idea of a pre-configured empty channel 
is hypothetical. Conformational flexibility induced by 
strong binding of the smaller Na+ in the selectivity filter 
is central to the mechanism of Na+ discrimination. 

This work supports the view that selectivity can be ad- 
dressed by analysis of local interactions involving a single- 
ion binding site j4i]. The assessment of rates, though, 
likely requires an account of interactions between ions 
occupying different channel binding sites [19|, but fluc- 
tuations of the sort identified here probably play a role 
too. 

The identification of channel-reorganization (induced 
fit) fluctuations is suggestive of solvent reorganization in 
chemical dynamics. The fact that a simple distribution, 
gaussian with slight positive skewness as in Fig. |31 works 
satisfactorily will have important conceptual and practi- 
cal consequences for modeling aqueous electrolyte solu- 
tions more broadly. 
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Figures 




FIG. 1: Depiction of the channel filter. Only two of the four 
chains are shown for clarity. Heavy atoms in the blue rib- 
bon segment experience a harmonic mean field, i.e., they are 
"restrained" with spring of force constant fc = 10 kcal/mol- 
A^. The segment colored yellow is the TTVGYG sequence 
of amino acids comprising the selectivity filter. The S2 site 
defined by VG is shown with the bound ion. 
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FIG. 2: An example how coupling an ion and the filter 
might enhance fluctuations. Here a coordinate in a parabolic 
potential-energy (— ■ — ) well also experiences a coulombic at- 
traction ( ) displacing the minimum energy position left- 
ward (unbroken line) . The fluctuations observed for the fully- 
coupled system, corresponding to the potential energy func- 
tion given by the solid curve, characterizes the available states 
on the harmonic potential function in the neighborhood of the 
physical binding energies. The curvature at the minimum is 
reduced, fluctuations are enhanced, and the observed fluctu- 
ations raise the free energy above the mean binding energy, 
according to Eq.|3 




FIG. 3: Normalized distribution of carbonyl oxygens radially 
from the ion, for the k = (hybrid) case described in the text 
for the wild type (bold lines) and the (GG) mutant (lighter 
lines). Notice that the longest ion-oxygen distance is about 
the same in each case, but the shortest ion-oxygen distances, 
and the most probable ones, are shorter in the case of Na+. 
This suggests that Na+ isn't merely delocalized, but constricts 
the filter relative to the K"*" case. 
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FIG. 4: Normalized radial distribution of 'body-diagonal' OO 
pairs, indicated by the arrow, for the selectivity filter in the 
cases of Na"*" and K"*", indicating relative constriction of the 
filter by Na"*". Bold lines: wild type. Lighter lines: (GG) 
mutant. Inset: atoms of the selectivity filter, showing an 
example of the 00 body-diagonal. Carbonyl oxygen atoms 
are in red, carbonyl carbons are green, and the ion is orange. 
Notice that the width of the distribution for the Na^ case is 
slightly the larger, and that relative widths, width/rynsx, of 
these distributions are roughly 10%. 
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FIG. 5: Normalized distribution of ion carbonyl binding en- 
ergies for the k — case for the wild-type. The circles (K^) 
and triangles (Na^) are the raw data. The binding energy 
abscissa is measured relative the mean binding energy, (e), 
for the Na"*" case. Thus, the mean binding energy for K"*" is 
about 22-25 kcal/mol higher than for Na"*". Gaussian fits to 
the data are smooth curves, and the additional dotted line on 
the left superposes the fit for the case on the Na"*" data. 
This shows that distribution is distinctly broader for the Na^ 
case. Similar results are obtained for the k = 10, 5.0, 2.5 cases. 
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FIG. 6: Free energy results for several values of the force 
constant k used in the harmonic mean field. The lower gray 
band is A/x'"^ (filter) of Eq. Q. The upper gray band is Afi'"^, 
the left-side of Eq. The lower boundary of each grey 

band is the result of the gaussian model Eq. |3] and the upper 
boundary was obtained on the basis of coupling-parameter 
integration. Note that the mean interaction energy differences 
(circles) are negative, signifying the better interaction of Na^ 
with the channel. 



